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Two-magnon Raman scattering provides important information about electronic correlations in 
the insulating parent compounds of high-T c materials. Recent experiments have shown a strong 
dependence of the Raman signal in Bi g geometry on the frequency of the incoming photon. We 
present an analytical and numerical study of the Raman intensity in the resonant regime. It has been 
previously argued by one of us (A.Ch) and D. Frenkel that the most relevant contribution to the 
Raman vertex at resonance is given by the triple resonance diagram. We derive an expression for the 
Raman intensity in which we simultaneously include the enhancement due to the triple resonance 
and a final state interaction. We compute the two-magnon peak height (TMPH) as a function of 
incident frequency and find two maxima at Wres ~ 2A + 3 J and cj'S ~ 2A + 8 J. We argue that the 
high-frequency maximum is cut only by a quasiparticle damping, while the low-frequency maximum 
has a finite amplitude even in the absence of damping. We also obtain an evolution of the Raman 
profile from an asymmetric form around w« s to a symmetric form around uJrea- We further show 
that the TMPH depends on the fermionic quasiparticle damping, the next-nearest neighbor hopping 
term t' and the corrections to the interaction vertex between light and the fermionic current. We 
discuss our results in the context of recent experiments by Blumberg et al. on Sr^CuOiCli and 
Y Ba2CusOe.i and Riibhausen et al. on PrBa2Cu^Oj and show that the triple resonance theory 
yields a qualitative and to some extent also quantitative understanding of the experimental data. 



I. INTRODUCTION 

In recent years a lot of efforts have been undertaken to 
understand the pairing mechanism in high-T c supercon- 
ductors ■ Some of the existing theories consider an 
effective electron-electron interaction mediated by spin 
fluctuations as the source of the pairing mechanism ||,|| ■ 
In the parent compounds of the high-T c materials the 
strong magnetic correlations lead to the occurrence of 
antiferromagnetism. Two-magnon Raman scattering is a 
valuable tool in probing antiferromagnetism and can thus 
provide important insight into the nature of the pairing 
correlations ||j9). 

The two-magnon Raman scattering cross section (Ra- 
man intensity) is proportional to the Golden Rule tran- 
sition rate JlCj ] 
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where i and / are the initial and final states of the sys- 
tem, 6ij are the corresponding energies and (e$ — ef) is 
the total energy of the two magnons in the final state. 
Mr =< e*f\Mn\ei > is the Raman matrix element (Ra- 
man vertex), ej and e/ are the polarization unit vectors 
of the incident and outgoing photons, and the summation 
runs over all possible initial and final electronic states. 

Graphically, the Raman intensity is given by the dia- 
gram shown in Fig. |l|a, where the intermediate magnons 
(wavy lines) are on the mass shell. The dashed lines in 
this diagram describe the incident and outgoing photons 



and the shaded circles represent the full Raman vertices, 
which include all effects of the final state magnon-magnon 
interaction ||llf . In conventional Raman experiments, one 
measures the Raman intensity, R, as a function of trans- 
ferred photon frequency Alu = uii — tOf where u>i and u) t 
are the incident and outgoing photon frequencies, respec- 
tively. The fingerprint of antiferromagnetism in these ex- 
periments is the presence of a two-magnon peak in R(Auj) 
|l2|. In the insulating parent compounds of the high-T c 
materials, this peak occurs at a transferred frequency of 
about 3000cm" 1 . The two -magnon peak has not only 
been observed in the insulating compounds, but also in 
electron 13 and hole doped materials |l4[ . 

Theoretically, most of the analytical and numer- 
ical studies of two-magnon Raman scattering were 
performed within the conventional phenomenological 
Loudon-Fleury (LF) theory |l5|| which assumes that the 
matrix element Mr for the interaction between photons 
and magnons is frequency independent. This implies that 
the theory neglects the internal structure which the ma- 
trix element possesses since the spin-photon interaction is 
actually mediated by fermions: the incident photon cre- 
ates a particle-hole pair which emits two spin excitations 
and annihilates into an outgoing photon (see e.g., Fig. |^). 
Despite this weakness, the LF theory was originally con- 
sidered as a suitable theory for Raman scattering in the 
parent high-T c materials because it predicts that the two- 
magnon profile should have a peak at a transferred fre- 
quency of about 2. 8 J where J is the magnetic exchange 
interaction. A comparison with the data |T^,|l6|,[l7j then 
yields J = 0.12eV which is fully consistent with the value 
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FIG. 1. (a) The Raman cross section is proportional to 
this diagram where the intermediate magnons (wavy lines) 
are on the mass shell. The dashed lines present incoming and 
outgoing photons and the filled circles are the full Raman 
vertices M R ot . (b) The full Raman vertices include all effects 
of multiple magnon-magnon scattering. The open circles are 
the bare Raman vertices Ma and the open squares describe 
the magnon-magnon scattering vertex V(k,q). 



for the in-plane exchange interaction extracted from neu- 
tron scattering (T|| and NMR data ]l9) . 

Recent experiments on single-layer SriCuOiCl2 
and double-layer YBa^Cu-iO^.x [^0| as well as on 
PrBaiCu-j.O'j plJl , however, presented some qualitative 
features of the Raman signal which cannot be explained 
within the framework of the LF theory. In these ex- 
periments, the Raman intensity was measured both as 
a function of transferred frequency at a given incident 
frequency tOi (the two-magnon profile), and as a func- 
tion of tOi at a fixed transferred frequency Alo ~ 2.8 J at 
which the two-magnon profile exhibits a maximum. In 
the latter case one in fact measures the variation of the 
two-magnon peak height (TMPH) with u>i. The experi- 
mental features which are in disagreement with the LF 
theory include: 

1) A strong dependence of the TMPH on a;, with two 
distinct maxima at tOi = coiW ~ 2 A + 3 J and at u>i = 

(2) 

cures ~ 2A + 8J, where 2 A ~ 1.7 'eV is the charge transfer 
gap Despite quantitative differences between vari- 

ous compounds, the second maximum in all compounds 
is always stronger than the first one. The LF theory, on 
the contrary, predicts that the intensity should only un- 
dergo a weak (logarithmical) enhancement at u>i = 2A 
and uji — 2A + 2.8J (ingoing and outgoing resonances). 
No enhancement, however, has been experimentally ob- 
served at u>i = 2A. 

2) The shape of the two-magnon profile is asymmetric 
and possesses a shoulder-like feature for transferred fre- 
quencies above the two-magnon peak, i.e., for Alu > 2.8 J. 
This feature has been observed around the first resonance 
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it practically disappears when the frequency of 



FIG. 2. The triple resonance diagram which yields the 
dominant contribution to the Raman intensity in the reso- 
nant regime. Solid and dashed lines represent fermions from 
conduction and valence bands, respectively. Notice that this 
diagram contains intraband scattering at the fermion-magnon 
vertices. 



the incident photon approaches the second resonance at 

LU res . 

Motivated by these findings, several groups studied 
two-magnon Raman scattering beyond the LF approx- 
imation p3|-p5[. It has been shown that the validity of 
the LF theory is restricted to the nonresonant regime, 
when the frequency of the incident light is much smaller 
than the charge-transfer gap 2 A |2aJ2l]. Most of the 
experiments, however, are performed with photon fre- 
quencies slightly above the charge transfer gap. In this 
resonant regime the internal structure of the Raman ma- 
trix clement cannot be neglected. Chubukov and Frenkel 
(hereafter referred to as CF) developed a diagrammatic 
approach to Raman scattering in the framework of the 
large-C/ spin-density-wave (SDW) approach to the Hub- 
bard model at half-filling 24 1. They identified those di- 
agrams which reproduce the LF vertex, and in addition 
identified a new diagram which is not included in the LF 
theory but yields the dominant contribution to the scat- 
tering process in the resonant regime. This new diagram 
has the largest amplitude when \luij — 2A| = O(J), and 
in addition, diverges in the absence of a fermionic damp- 
ing in a small region (nearly a single critical line) in the 
(wi, Alu) plane where all three terms in the denominator 
vanish simultaneously (see Fig. ||). Due to this property, 
the new diagram identified by CF is called the triple res- 
onance diagram. The inclusion of a fermionic damping 
eliminates the divergence, but the Raman matrix element 
remains strongly peaked along the critical line. Since 
the computation of the full Raman intensity with the 
triple diagram for the Raman vertex is rather involved, 
CF used a semi-phenomenological approach to analyze 
the dependence of the TMPH on cjj. They considered 
the final state magnon-magnon interaction and the triple 
resonance enhancement separately, and conjectured that 
the experimentally observed two maxima in the TMPH 
occur at LUi for which the Raman vertex resonates at the 
same transferred frequency Aw = 2.8J at which the two- 
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Cfli - 2A in units of J 

FIG. 3. The shaded area represent the region in the 
(wi,Aa>) plane in which the triple resonance occurs. The 
horizontal line corresponds to Aui = 2.8J at which the 
two-magnon profile has a peak. 



magnon profile has a peak. By analyzing where the reso- 
nant line for the Raman vertex crosses Aw = 2.8 J, they 
obtained two resonance frequencies, LuiW ~ 2.9 J + 2A 
and ~ 7.9 J + 2 A which both agree with the experi- 
mental data. 

The analysis in Ref. p4| , however, left several issues 
open. First, the validity of the semi-phenomenological 
approach needs to be verified. Second, the quantita- 
tive behavior of the TMPH as a function of iOi and, in 
particular, the form of the peaks at oj^es and w„ s and 
their relative amplitudes have not been studied yet. CF 
merely conjectured, without performing explicit calcula- 
tions, that the resonance at Ljies should be weaker than 



the one at w^el because near ct^ei) there exists a strong 
restriction on the possible directions of the magnon mo- 
menta which satisfy the resonance condition at a given 
magnon energy. No such restriction exists near ajf-es- 
This conjecture also has to be verified by explicit cal- 
culations. Third, the anisotropy of the two-magnon pro- 
file and its evolution with varying incident frequency has 
not been studied. Forth, the calculations in CF were per- 
formed in the framework of a mean-field, large U , spin- 
density-wave (SWD) approach to the Hubbard model 
with only nearest-neighbor hopping. This theory, how- 
ever, possesses the weakness that it predicts that the 
maximum of the valence band is degenerate along the 
boundary of the magnetic Brillouin zone. Meanwhile, 
experiments on SriGu02Cli have demonstrated that 
the valence fermions possess a strong dispersion along 
the magnetic Brillouin zone boundary with maxima at 
(7r/2, 7r/2) and symmetry related points |2^,|2^]. This dis- 
persion can easily be reproduced in the SDW formalism 
if one includes a next-nearest-neighbor hopping, t' . This, 
however, changes the energy denominator in the triple 
resonance diagram, and one therefore has to reexamine 
the conclusions of CF by performing their calculations 



for the (t — t' — U) model. The inclusion of t' is par- 
ticularly relevant for computations near tuiW since the 
dominant contribution to the Raman vertex in this fre- 
quency range comes from fermions near the top of the 
valence band whose degeneracy is lifted by t' . 

The goal of the present paper is to address the above is- 
sues. We will compute below the Raman intensity includ- 
ing both a final state interaction and the enhancement 
of the Raman vertex due to the triple resonance. We 
will study the Raman profile and the TMPH numerically 
and analytically and demonstrate that the two peaks in 
the TMPH survive the effects of the magnon-magnon in- 
teraction. We will analyze the relative amplitude of the 
TMPH near u>ies and Wrel and show that although the di- 
vergent piece near uii-es is much weaker than the one near 
LOrJs, the nondivergent term is much larger near u>res- As 
a result, the relative amplitude of the two peaks in the 
TMPH turns out to be strongly dependent on the quasi- 
particle damping which cuts the divergent part but does 
not affect the subleading term substantially. We will also 
show that both peaks in the TMPH are anisotropic - the 
intensity drops much faster on the high-frequency side 
of each of the peaks. Further, we will study the effects 
on the TMPH of a next-nearest neighbor hopping t' and 
vertex corrections to the interaction between light and 
fermionic quasiparticles. 

We will also study how the two-magnon lineshape 
evolves with the incident frequency, and show that it 
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changes from an asymmetric form for u>i 
symmetric form around u>i « w^el- 

We will compare our results with the experimental 
data on on Sr2Cu02Cl2 and Y BCI2CU3O6.1 by Blum- 
berg et al. and on PrBa2CuzOi by Riibhausen et al. 
and demonstrate that all the features in the two-magnon 
profile and the TMPH observed in Raman experiments 
can be qualitatively described by the triple resonance di- 
agram. At the same time, we will see that quantitative 
agreement with the data is not always perfect. 

The paper is organized as follows. In Sec. |n], we present 
the formalism and the expressions for the Raman vertex 
both in the LF approximation and near the resonance. In 
Sec. ra we present our analytical results for the vertex 
near Uri'P and demonstrate that the divergence is much 
stronger near the upper resonance frequency. We also 
discuss in this section how the inclusion of a next- nearest- 
neighbor hopping t' affects the resonance behavior of the 
Raman vertex. In Sec. IV we present our numerical re- 



sults for (a) the Raman line shape for different incident 
frequencies (Sec. IV A| ) and (b) for the TMPH as a 
function of ojj (Sec. IV B|). We then discuss its depen- 



dence onthe fermio nic d amping and the inclusion of t' 
(Sec. [V C ). In Sec. |IVD| we consider vertex corrections 
to the interaction between light and the fermionic cur- 
rent. Finally, in Secjv] we compare our results with the 
experimental data. 
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II. THE FORMALISM 



We can now decompose the interaction vertex into 



The two-magnon Raman scattering cross section is 
given by Eq.(|l]). In this paper, we will focus on the x'y' 
scattering geometry where most experiments have been 
performed. In this geometry, the polarization unit vec- 
tors of the incident and outgoing light are both real (lin- 
early polarized light), perpendicular to each other and di- 
rected at 45° to the crystallographic directions, i.e., ii = 
x 1 = (x + y)/ V2, if = y' = (x — y)/ \f2. Other scattering 
geometries for linearly polarized light are x'x' where ii = 
if = (x + y)/ V2, xy where ii = x,if = y and xx where 
ii = x, if = x. For circularly polarized light, the scat- 
tering geometries are LL where ii = if = (x + iy)/ V2, 
and LR where ii — (x + iy) / V2, if — (x — iy)j\f2. 

The x'y' geometry is often referred to as the B\ g ge- 
ometry. Strictly speaking the matrix element in the 



x'y' scattering channel includes B\ g 

r-Bi 



and A2 g compo- 



nents with M 
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x\M R \x > - < y\M R \y > and 

< x\Mjt\y > — < y\Mn\x >. In practice, how- 
ever, the Ai g component is always negligible, and we will 
therefore not distinguish between x'y' and pure B\ g scat- 
tering. 

As we already discussed in the introduction, we will 
consider Raman scattering in the framework of the large- 
ly SDW formalism for the one-band Hubbard model. In 
this approach, one introduces a long-range antiferromag- 
netic order and decouples the electronic dispersion into 
two subbands of valence and conduction fermions. 

The diagrammatic approach to the Raman scattering 
in the SDW formalism was developed by CF. An exam- 
ple for the diagrams which contribute to the bare Ra- 
man vertex is shown in Fig. |[ This diagram contains 
two types of vertices: one for the interaction between 
fermions and light, and one for the interaction between 
fermions and magnons. The interaction with light ap- 
pears in the SDW theory as a result of the modulation 
of the hopping matrix element by the vector potential of 
the electromagnetic field. The spin-fermion vertices can 
be straightforwardly obtained from the full expression of 
the spin susceptibility which in the SDW theory is given 
by the RPA series of bubble diagrams. In each of theses 
bubbles one fermion is from the valence band and the 
other is from conduction band. 

To obtain the Raman intensity, we need to know 
the full Raman vertex which includes a whole series 
of magnon-magnon interaction events. It has already 
been emphasized several times in the literature that the 
dominant contribution to the magnon-magnon scattering 
comes from the region near the magnetic Brillouin zone 
boundary where the antiferromagnetic magnons behave 
almost as free particles |11 2q . In this situation, the only 
relevant interaction term has two creation and two anni- 
hilation magnon operators: 
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Vk-q = VkVl + VkVl + VkVl + VkVl . (3) 

where the different symmetry factors are given by 

Vk = ~(cosfc x + cos/c^); v k = i(cosfcj; - cos k y ); 

D k = ^(sinfca; +sin/c y ); u k = ^(sink x -sink y ). (4) 

Before we discuss our calculations for the full Raman in- 
tensity at resonance, we briefly review the calculation of 
the Raman intensity in the nonresonant regime when the 
LF theory is valid. In the LF theory, the bare Raman 
vertex (open circle in Fig. |l|b) is assumed to be indepen- 
dent of the photon frequencies while its dependence on 
the magnon momentum q has the form |23 24] 



Mr = A v q (e lx e* fx + e iy e fy ) 

-(e ix e*f x cosq x + ei y e fy cose 



where A is a constant. In the diagrammatic approach, 
the LF vertex is obtained by collecting the diagrams with 
interband scattering at the magnon-fermion vertices. At 
photon frequencies small compared to the SDW gap, 
these diagrams have the largest overall factor. One can 
easily check that Mr is finite only in the B\ g scattering 
geometry and for LR polarized light. In both cases we ob- 
tain Mr — ~AD q . For this particular form of Mr, only 
the second term in Eq.(|3]) contributes to the magnon- 
magnon scattering process. With this simplification, the 
summation of the ladder series for the full Raman in- 
tensity can be reduced to solving an algebraic equation. 
Doing this we obtain for the full Raman intensity in the 
B\g channel 



R(uj) oc Im 



1 + //4S 

where S is the value of the spin, and 

AJStt-^ (cosq x — cosq y ) 2 
Auj — 2uj q + i5 



1J J X ^ | 



(5) 



(6) 
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with Ao; = u>i — u!f and uj q — AJSyl — is q is the magnon 

dispersion. Eqs.(|5|) and (^) yield a two-magnon peak at 
Aui = 2.8J (for S = 1/2), but R(ui) clearly contains no 
dependence on the incident photon frequency oji |]IT| , |l5| . 
As mentioned earlier, the LF theory is only valid for small 
u>i. When oji is comparable to the gap between the con- 
duction and valence bands (which in the cuprates is the 
charge transfer gap) , it turns out that diagrams with in- 
traband scattering at the fermion-magnon vertices (in 
contrast to interband scattering in the LF diagrams) be- 
come dominant. The most relevant of these diagrams is 
shown in Fig. ^. This diagram is called the triple res- 
onance diagram because it contains three terms in the 
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FIG. 4. The effective Raman vertex incfudes a single 

magnon-magnon scattering event. 



denominator which can all vanish simultaneously if we 
adjust the incident and final photon frequencies. The 
analytical expression for this diagram is given by 
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and the prime indicates summation over the magnetic 
Brillouin zone. The iT term represents the fermionic 
quasi-particle damping which we assume for simplicity 
to be independent of momentum. The actual damping, 
indeed, should have some momentum dependence, par- 
ticularly near the top of the valence band [p9| |32|| . Out 
of the three terms in the numerator, the first two are 
the vertex functions for the interaction between light 
and fermions, while the third term is the product of 
the two vertices for the interaction between fermions 
and magnons. As follows from Eq.(Q), the resonant Ra- 
man vertex Mr depends on the magnon momentum, but 
also, via the denominator, on the incident and outgoing 
photon frequencies and on the magnon frequency. The 
frequency dependence of Mr can only be eliminated in 
the artificial limit when the fermionic damping is very 
large and overshadows all other terms in the denomi- 
nator. For small and moderate damping, the frequency 
dependent and also momentum dependent terms in the 
denominator cannot be ignored, and this makes the com- 
putation of the full Raman intensity rather involved. We 
found, however, that the diagram which is most difficult 
to compute is the one without a final state interaction. 
At the same time, the series of diagrams with at least 
two scattering events can easily be summed up because 
the Raman vertex renormalizcd by the inclusion of just 



a single magnon-magnon scattering event no longer de- 
pends on the magnon frequency while its dependence on 
the external magnon momentum I reduces to a simple 
i>i form for B\ g scattering. This effective Raman vertex, 
M R ff = vi M R ff is shown in Fig. |. 

We remind that the experimentally measured Ra- 
man profile for any incident photon frequency contains 
a prominent two-magnon peak which is solely due to 
magnon-magnon scattering. In this situation, the dia- 
grams without and with a single magnon-magnon scat- 
tering are most likely to be less relevant than the dia- 
grams with multiple scattering events. For our analytical 
considerations, we neglect the diagrams with no or only 
one magnon-magnon scattering event. In this approxi- 
mation, we can formally rewrite the full Raman intensity 
in the same form as for the LF theory: 



R(u>i,u>f) oc Im< 
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where / is the same as in Eq. (|), and M R ff is obtained 
by substituting Eq.(^) into the diagram in Fig. | and 
performing the integration over the intermediate magnon 
frequency and momenta. We then obtain 
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(^Auj — 2uj q + iS^j (uii — uj q — E k — E k - q + iT 

It is essential however that M R ^ still possesses a complex 
dependence on the external incident and outgoing photon 
frequencies. This in turn implies that the full intensity 
R(LUi,ujf) is a function of both frequencies rather than of 
LUi — ujf as in the LF theory. A very similar approach was 
used by Schonfcld et al. [p5| . 

In our numerical calculations of the full Raman inten- 
sity we considered all diagrams, i.e., diagrams with zero, 
one or multiple magnon-magnon scattering events. The 
details of this computation are presented in Appendix 
We found a good qualitative agreement between our 
numerical and analytical results and consider this as a 
partial justification for the omission of the lowest order 
diagrams in our analytical considerations. 

We now proceed with the discussion of our analytical 
results, and then present our numerical data. 



III. ANALYTICAL RESULTS 

In this section, we present the results of our calcula- 
tions of the Raman intensity near the two resonant fre- 
quencies, UJres and We first consider the case t' = 
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conduction band 
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FIG. 5. Quasiparticles at the bottom of the valence band 

(2) 

contribute to the Raman vertex at Ures, whereas quasiparti- 
cles at the top of the valence band contribute to tJrei- The 
dashed and wavy lines represent the excited quasiparticles 
and the emitted magnons, respectively. 



and then discuss how the intensity changes if we break 
the particle-hole symmetry by including a hopping term 
between next-nearest neighbors. Our point of departure 
is the approximate expression for the Raman intensity, 
Eq.(||), in terms of the effective Raman vertex, M^. 
We first study the form of M e ^ near w™' and then dis- 



cuss the form of the vertex near u> 



(i) 



(2) 



A. Resonance at u, 

As we discussed in the introduction, the upper reso- 

(2) 

nance frequency, Wres, is close to the maximum possible 
incident frequency 2A + 8 J for which the resonance con- 
dition can still be satisfied (we assume here that A is 
large, i.e., E/. = A + J(cos k x + cos k y ) 2 ). The resonance 
at the highest incident frequency corresponds to a pro- 
cess in which light causes a transition of a quasiparticle 
from the bottom of the valence band to the top of the 
conduction band (see Fig. ||). Since the bottom of the 
valence band is at k = (0, 0) and is not degenerate at the 
mean-field level, it is reasonable to assume that the cor- 
rections to the mean-field fcrmionic dispersion near the 
bottom of the band will only give rise to a finite lifetime, 
but will not introduce any new qualitative features into 
the fcrmionic spectrum. We therefore restrict our con- 
siderations to the mean-field form of the fermionic dis- 
persion and model the fluctuation effects by introducing 
a quasiparticle damping T . 

To simplify the presentation, we first neglect the nu- 
merator in Eq.(^|) and focus on the resonance behavior 
of the effective vertex due to the vanishing denomina- 
tor. Near the bottom of the band, we can expand the 
fermionic dispersion and obtain Ek = A + 4 J — 2Jk 2 . 



We will assume that near the resonance the bosonic mo- 
mentum q is also small such that the expansion near the 
bottom of the band holds for both Ek and Ek+ q . This 
last assumption will be verified after we perform our cal- 
culations. Further, for small q we can expand uj q as uj q — 
J\^2q. Substituting the expanded forms of Ek, Ek+ q and 
uj q into Eq.(^|), neglecting the numerator, and introduc- 
ing the dimensionless variables A^ = (u>i — 2A)/2J and 
A A = Aw/2 J we obtain 



M e J S oc 



d 2 k d 2 q 



{AX-qV2 + iT) 
1 



(Ai - 4 + 2k 2 + iT)(\i - 4 - AA + 2{k - q) 2 + iT) 
1 

(2A, - 8 - AA + 2k 2 + 2(k - q) 2 + (AA - qV2) + iT) ' 

(10) 

A similar expansion has been performed by CF. They 
however considered only the bare Raman vertex (the 
one without a final state interaction) in which case the 
Golden Rule requires the magnons to be on the mass 
shell, i.e., q = qa = AX/y/2. For on-shell magnons, the 
integration over d 2 q just yields a factor 2ir 2 iqo. Further- 
more, for q = qo, the last term in the denominator in 
Eq.(|l0[) is the sum of the second and the third term, 
so one can tune k and the angle between k and q such 
that all three terms in the denominator of ([ll]) vanish 
simultaneously. Expanding near this point, CF obtained 
M fl cx(A,-A re;5 )- 3 / 2 . 

We now show that the fully renormalized M^' pos- 
sesses the same functional form as the bare vertex ob- 
tained by CF since the integration over q by itself is 
confined to the vicinity of q$. To demonstrate this, we 
expand near the point where the second and the third 
term vanish and then integrate over k and q. The 
conditions that the two denominators vanish simulta- 
neously at a given q are 2k 2 — 2k 2 = 4 — A^ and 
cos</> = cos (fo = (2q 2 — AA)/(4fc g)- We will see below 
that in order to obtain the most singular contribution to 
one has to expand around ko and 4>q to linear or- 
der in x = k — kg and to quadratic order in y = (f> — </V 
Performing this expansion, we obtain 



where 
G(q) 



M e R ff oc 



dx 



d 2 q G(q) 



AX-qV2 + iT 



1 



(11) 



dy 



(x + iT){k x + (AA - g\/2)/4 + iT) 
1 



C\ + fc o gsin0o2/ + k qy 2 /2 + iT 



C 2 + k qsm4> y + k qy 2 /2 + iT 



(12) 



G 



and G\ and C2 are given by 

C x = x(k - {2q 2 - AX)/4k ) ; 
C* 2 = d + k x + (AA - gV2)/4 . 

The integration over y can be done explicitly and yields 
G(q) = 2ir ' 



(x + iT)(k x + (AA - <7\/2)/4 + iT) 
1 



L \l 2fc gCi - (k q) 2 sin 2 O + «r 
1 

2k qC2 - (koq) 2 sin 2 O + iT 



(13) 



We immediately see that if sin <p is finite, one can expand 
the square root in the integrand in Eq.(j^) in C1.2 and 
find after simple manipulations that the Raman vertex 
is free from singularities. To obtain a divergent, reso- 
nant contribution to the vertex, we therefore have to set 
sin <pQ = 0. For any given q, this singles out a line in the 
(Aj, AA) plane with 



A i = A(, 2 )( g ,AA) 



1 



AA 



(14) 



Furthermore, we also have to require that the poles and 
branch cuts in the integrand in Eq.(|l3|) be in different 
half-planes, since otherwise, the integral over x just van- 
ishes. This requirement is satisfied if C1.2, which both are 

linear functions of x, have negative derivatives dCi^/dx. 

(2) 

Near the resonance line A; = Ares, we have 



Ci 
C 2 



(AA) 2 



x 4 
"Ifo ~ Aq 2 
x {2q 2 - AA)AA AA - y/2q 



2q 2 



We see that both derivatives are negative if 2q 2 > AA. 
If, as we assume, q is close to q = AA/V2, then the 
derivatives are negative provided A A > 1. For trans- 
ferred frequencies near the two-magnon peak, we have 
A A ~ 1.4, i.e., the above condition is satisfied. Perform- 
ing then the integration over x, we obtain keeping sin <pQ 
small but finite 



G(q) = 



16tt 2 



1 



AA - qy/2 k q 



A\-qV2 
2k q 



iT 



sin 0o 



A\-qV2 I 
2k q I 



I sin 00 I 



(15) 



We now express sin 2 (j)o in terms of deviation from the 
critical line as 



sin 2 0o = S - (AA - qV2) 



AA + 1 
AA(AA - 1) 



(16) 



where we introduced <5 = A, — A 



(2) 



We now substitute 



( |l5| ) and ( [l6| ) into Eq.(|Tl]) and perform the integration 
over q. Introducing z — (AA — qV%) /$ we finally obtain 



M e J f oc 



Z(AA) 1 
AA <5 3 /2 ' 



(17) 



where 



Z(AA) =2(1 + AA) f 

J — r 



dz 



-00 (z-X) 2 

1 



vT 



1 



z-iV 



2(1+AA) 



is a very smooth (nearly a constant) function of the mo- 
mentum transfer. Clearly, the integral over z comes from 
2 = 0(1) which in turn implies that the actual integral 
over q is confined to the region AA - qV2 = 0(6). The 
fact that q is confined to the vicinity of its on-shell value 
implies that one can set q = AX/^/2 in Eq.(p^) which 
yields 



A r 2 2(AA) 



(AA-1) 2 



(18) 



For AA = 1.4, which, we remind, corresponds to the po- 

(2) 

sition of the two magnon peak, we obtain A res = 3.93. 
Also, for this AA we have q w 1. This q is not small 
enough to fully justify our expansion in the magnon mo- 
mentum. However, for q x = q v , the approximation of cu q 
by a linear term is incorrect only by 8%. 

We see that Mj' diverges as 5~ 3 / 2 when the incident 

(2) 

frequency approaches the resonance value A res . The same 
result has been obtained by CF who neglected the final 
state interaction. In this respect, our result shows that 
the final state interaction does not destroy the triple reso- 
nance which occurs along the same line as in the absence 
of the magnon-magnon scattering. 

We now consider the effect of the numerator in M|j f. 

(2) 

In principle, the numerator is finite at \ = A res so that 
the <5 -3 / 2 behavior should survive close to the resonance 
line. In practice, however, for transferred frequencies 
around AA = 1.4, the resonance value of Ai is close to the 
maximum possible resonance value of 4. At this maxi- 
mum value of Ai the matrix element for the interaction 
between light and the fermionic current vanishes and the 
numerator turns to zero. A simple analysis similar to the 
one performed by CF shows that the numerator vanishes 
as 4 — A,- for A 



4. The full M e ^ s then behaves as 



Mr* = 



A(AX) 



A, 



|A 4 -A$4i|3/2 



(19) 



Restoring all numerical factors in the expression for the 
Raman vertex, we find that A(A\) has the form 



A(AX) 



J ( t 

7T \2J 



(1 + AA) Z(AX) 



(20) 



7 



(2) 

At some distance away from Ares , the difference between 
4 — Xi and Xr es — \ is irrelevant and M.^ can be ap- 
proximated by 



M 



eff 



A(AX) 



IA (2) 

| s^res 



(21) 



This is the final expression for the effective Raman ver- 
tex. Substituting this expression into Eq. (||) , one obtains 
for the Raman intensity 



R(X h AX) oc 



\A(AX)f 



I 



(2) 
res 



A,: | 1 + //4S 



(22) 



where / is given in Eq.(|). This result for R is in fact 
very similar to what CF conjectured: the full intensity is 
a product of two terms, one, 1/(1 + I/AS), depends on 
AA and has the same form as in the LF theory, while the 
other, |Aj — X^J S in essence reflects the enhancement of 
the bare Raman vertex. This last term depends on Xi but 

(2) 

also on the transferred frequency through A res . There is 
however an extra feature: the interplay between triple 
resonance and magnon-magnon interaction gives rise to 
an extra dependence of the intensity on AA, through the 
factor A. However, we found that A is a smooth function 
of momentum transfer, so this extra dependence is not 

(2) 

that relevant. This is particularly true near A res where 
the two-magnon peak is rather narrow so one probes A 
only in a narrow range of AA around 1.4. We already 
mentioned in the introduction that the inverse linear de- 

(2) 

pendence of the TMPH near u>res was first experimentally 
observed by Blumberg et al. for Y Ba^Cu-iO^.x [^0| and 
later by Rubhausen et al. for PrBa^Cu-iOi pip . This is 
fully consistent with our result. We, however, performed 
our calculations only in the vicinity of the resonance and 
therefore cannot provide a theoretical estimate for the 
width of the region where the inverse linear behavior 
holds. The data for Sr2CuC>2Cli are a bit less conclusive 
because in this material, w res is larger than the largest 
experimentally accessible Wj, and one cannot unambigu- 
ously conclude from the data that the Raman intensity 
follows an inverse linear behavior in a wide range of fre- 
quencies. 

Besides the behavior of the TMPH near A r |i, Eq.(§|) 
also describes the form of the two-magnon profile at a 
given Xi. The conventional factor 1/(1 + 1 /AS) produces 
a symmetric peak in the intensity at AA = 1.4. The other 
two factors which contribute to the peak lincshape are the 

(2) 

dependencies on A A in the overall factor A and in A res . 
Near the resonance value A^ = 3.93, we found that the 
two extra contributions to the lineshape cancel each other 
and the resulting two-magnon profile is chiefly given by 
1/(1 + I /AS) and therefore symmetric. For incident fre- 
quencies smaller than Xi = 3.93, the triple resonance in 
the Raman vertex occurs at transferred frequencies larger 
than A A =1.4. This obviously gives rise to an asymme- 
try of the two-magnon lineshape with a larger intensity 



at higher frequencies. The evolution from an asymmetric 
to a symmetric form of the two magnon profile when Xi 
approaches the resonance value from below is consistent 
with the experimental data. We will also demonstrate 
this effect when we discuss our numerical results. 



B. Resonance near u) r Vs 

The triple resonance theory of CF predicts that the 
TMPH measured as a function of LOi = 2A+2 JXi exhibits 
a second maximum at a relatively small frequency u) r ]) a = 
2 A + 3 J (see Fig. |J). For this low frequency resonance 
CF found a much smaller range of magnon momenta for 
which the resonance conditions are satisfied and therefore 
concluded that this resonance should be weaker than the 

(2) 

one at u! res ■ We now discuss this issue in more detail and 
will show that while the divergent term at LuiW is almost 
completely suppressed, the subleading terms are larger 
than at uji 2 e \. 

In the previous subsection, we expanded the fcrmionic 
dispersion around the bottom of the valence band since 

(2) 

LOres IS close to the maximum possible incident frequency 
for the triple resonance. Here, on the contrary, we will 
make use of the fact that the low-frequency resonance 
occurs at uj r \\ ~ 2A -1- 3 J which is not far from the min- 
imum incident frequency (= 2A + 2 J) for which a triple 
resonance is possible. Accordingly, we will study the be- 
havior of the Raman intensity by expanding the fcrmionic 
dispersion upto quadratic order around the top of the va- 
lence band. As in the previous subsection, we will also 
expand the magnon dispersion to linear order in the mo- 
menta. We will see that in this approximation the Raman 
vertex is free from actual divergencies. 

A peculiarity associated with the expansion of the 
fcrmionic dispersion near the top of the valence band 
is that the position of the band maximum is degenerate 
at the mean-field level. Numerous analytical and nu- 
merical calculations, however, have demonstrated that 
this degeneracy is an artifact of the mean-field treatment 
0,|9|j33|-|3T 



the actual fermionic dispersion possesses 
a maximum at (tt/2,tt/2) and symmetry related points. 
Recent photoemission experiments on Sr^CuO^C^ con- 
firmed this result and in addition have shown that the 
dispersion near the top of the valence band is nearly 
isotropic around (tt/2,tt/2) p6| , p7| ]. At the moment it is 
still a topic of controversy, whether one needs a substan- 
tially large next-nearest-neighbor hopping to explain an 
almost isotropic dispersion, or whether it is a property of 
the nearest-neighbor model in the strong coupling limit, 
as was first suggested by Laughlin |32| . The set of param- 
eters suitable to Sr^CuO^Cl^, namely J ft ~ 0.4, corre- 
sponds to an intermediate coupling regime in which case 
the next-nearest-neighbor exchange is probably needed 
to account for the isotropy of the spectrum |}(|. This 
second-neighbor hopping breaks the particle-hole sym- 
metry and causes substantial complications for the cal- 
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culation of the Raman vertex since Eq.([|) is no longer 
valid. At the same time, we know that the degeneracy 
along the magnetic Brillouin zone boundary is already 
lifted in the nearest-neighbor model due to self-energy 
corrections. One can therefore argue that the inclusion 
of t' only aids in fitting the ratio of the effective masses 
but does not introduce any new qualitative features. To 
avoid unnecessary complications, we will assume that the 
particle-hole symmetry is preserved, and that the experi- 
mentally measured nearly isotropic quasi-particle disper- 
sion around the top of the valence band is due to strong 
self-energy corrections. In other words, we will still use 
Eq.(^J) and will also assume that near the top of the va- 
lence band we can expand the fcrmionic dispersion as 
Ek = A + Jk 2 , where k measures the deviation from 
("7r/2, 7r/2). Substituting the expansion for Ek into Eq.(j^) 
and neglecting the numerator which is finite near the low- 
frequency resonance, we obtain 



M 



x 



d 2 k d 2 q 



eff 

11 ' 1 AX-qV2 + iT 
1 

x ■ 



(Aj - k 2 + iT)(Aj - AA - (k - q) 2 + iT) 
1 

(2A, - AA - k 2 - (k - q) 2 + (AA - qy/2) + iT) 



(23) 



To study whether the effective Raman vertex diverges at 
some particular A 2 ; , we perform the same analysis as in the 
previous section, i.e., expand near a particular fermionic 
momentum, ko where the second and third terms in the 
denominator vanish simultaneously at a given q. For ko 
and (f>o, which is the angle between ko and q we find 
feg = Ai and cosc^o = (q 2 + AA)/(2fco<z)- Expanding, as 
before, around ko and <po to linear order in x = k — ko 
and to quadratic order in y = <j) — <^0j we obtain 



M e J* oc I d 2 q 



G(q) 



where 

G(q) = I dx 



AA - qV2 + iT 



1 



(24) 



/ 



x / dy 



(x + iT) (k x + (AA - qV2)/4 + iT) 
1 



Ci + koqsm(f>oy - k Q qy 2 /2 + iT 
1 



C 2 + k qsin(poy - k qy 2 /2 + iT 
and C\ and G 2 are given by 

C x = x(k - (q 2 + AA)/2fc ) ; 
C 2 = Ci + k x + (AA - qV2)/4 



(25) 



The integration over y can again be performed explicitly, 
and we obtain 



G(q) = ~2m I dx 



1 



(x + iT)(k x + (AA - q\f2)/A + iT) 
1 



2fc qCi + (k q) 2 sin 2 (j} + iT 



2k qC 2 + (k q) 2 sin 2 4>o + iT . 



(26) 



This expression is similar to Eq.(|l3|) and we again find 
that if sin O is finite, one can expand the square root and 
obtain that G(q) is free from singularities. The expan- 
sion does not work, however, if sin^o = 0. In this case, a 
power counting argument indicates that the Raman ver- 
tex may diverge. The condition sin^o = singles out a 
line in the (Ai, A A) plane with 



A, 



= aW(< 7 ,aa) = <l 



1 



AA 

„2 



(27) 



The very existence of the critical line along which the 
Raman vertex diverges by the power counting argument, 
however, does not guarantee that the divergence is ac- 
tually genuine. Indeed, the arguments displayed in the 
previous subsection show that the vertex diverges only if 
the poles and the branch cuts in Eq.(p6|) are located in 
different half-planes. This is the case if the derivatives 
over x of C1.2 are negative. However, near Afel we have 



C 2 = Ci + xk 



A\-qV2 



Since (2 — AA) is always positive, the derivatives are 
clearly positive in which case the divergent contribution 
vanishes after the integration over x. We see therefore 
that the triple resonance does not yield a divergence in 
the Raman vertex at Ai = Aril, contrary to what we 
found near X^ s . 

We then studied the form of the Raman vertex in more 
detail and found that the absence of a divergence is a 
result of the restriction to a quadratic dispersion around 
the top of the band. Expanding further in k and redoing 
the calculations, we obtained a divergence in Mr result- 
ing from the integration over a small region of fermionic 
momenta. A similar result was also obtained by CF who 
used a somewhat different technique. However, the phase 
factor associated with the divergent contribution to Mr 
is very small, and the divergence is already eliminated by 
a small fermionic damping. 

So far, we have found that the Raman vertex exhibits 

a regular behavior around Aril. Experimentally, how- 
ever, the TMPH clearly displays a second maximum at 
Ai ~ 1.5. We will now show that this maximum can 
in fact be described within the triple resonance theory 
since the Raman vertex turns out to be very strongly 







enhanced near Aril- To demonstrate this, it is not suffi- 
cient to expand near particular values of fco an d 4>o, and 
we thus need to study the full form of the Raman ver- 

(2) 

tex. The full form of Mr near Ares was first obtained 
by CF and we now have to perform the same analysis 
near Afil. There is, however, a subtlety related to these 
calculations. CF had to assume that the magnons are 
on-shell, i.e., q = qo = AX/2, since a full analytical con- 
sideration is not possible without this last assumption. 
In the previous subsection, we found that the divergent 
piece of the vertex is intrinsically confined to a narrow 
range around q , and the results with and without the re- 
striction to only on-shell magnons are roughly the same. 
Near Afil, the situation is less rigorous since the diver- 
gent piece is absent. At the same time, it still looks 
reasonable to estimate the value of the vertex near Aril 
by just restricting with on-shell magnons. Doing this 
and following the computational steps outlined by CF, 
we obtain after some lengthy calculations 



M e R ff « S l 



|<J|3/2 



log- 



1 



'25 - iT 



'25 + iT 



log- 



1 



2\/2(5 
2+AA 



iT 



Do 



|«J|3/2 



log 



25 + iT 



1 - V25-iT 



log- 



2-AA 



2V2<5 
2+AA 

iT 



1 - P^r - iT 



2-AA 



iT 



(28) 



with Di = (2 + AA)/(\/2(AA) 2 ), D 2 = (2 - 

AA)/(\/2(AA) 2 ), and 5 = Afil — A, where Xres is given 
by Eq.(p7[) with q substituted by q : 



A«(AA) = i + ^ 



(29) 



Though both terms in Eq.(p8|) contain a term \5\ 3 / 2 , the 
combination of logarithms vanishes when <5 approaches 
zero. Moreover, expanding in <5, we find that the terms 
of 0(< 5 1 / 2 ) also cancel each other. Consequently, there 
is not even a weak triple resonance at Aril, and the Ra- 
man vertex turns out to be a regular function of Xi in 
the immediate vicinity of the would-be resonance line 
A r il ■ This indeed agrees with our expansion near /c and 
<f)Q. At the same time, it follows from Eq.(|2|) that at 

frequencies only slightly smaller than Aril, namely at 

(2 — AA)/2 < "v25, one of the logarithms contains an ex- 
tra in factor associated with the branch cut. Due to this 
extra factor, Mjj in fact scales as |<5| -3 / 2 , i.e., the vertex 
possesses the same functional dependence on the incident 
frequency as if the triple resonance at Afil were actually 
present. Near the two-magnon peak we have AA = 1.4. 
In this case the singular behavior actually starts very 
close to Aril, namely at Afil — A^ ~ 0.045. This singular 
behavior exists, with decreasing amplitude, upto 25 = 1, 
or Aril — Xi ~ 0.5, though at such high deviations from 



Aril the regular and singular parts of G(q) are of the 
same order. We see therefore that despite the absence 
of a true divergence at Afil, the TMPH still possesses 
a maximum very close to it. Moreover, for A A = 1.4, 
we have Afil = Afil ~ 1.45, i.e., u> re \ ~ 2A + 3J which 
is in good agreement with the experimentally observed 
location of the low-frequency peak in the TMPH. 
For experimental comparisons, it is essential that the en- 
hancement due to the branch cut in Mr is asymmetric - 



> A (1) 



This should 



it exists for A^ < Afil but not for A 
obviously yield an asymmetric form of the TMPH near 
Afil - the intensity should increase continuously as one 
approaches Afil from below and drop down rather fast 
when Xi exceeds Afil. In addition we should also obtain 
an asymmetry of the two-magnon lineshape at exactly 
Xi — Xres with a higher intensity at larger frequencies. 
Indeed, for AA > 1.4 we have Afil > A,-, and the Raman 
vertex is strongly enhanced. No such effect, however, ex- 
ists for AA < 1.4. Both above mentioned anisotropics 
are consistent with the experimental data. 

Finally, consider the numerator in the Raman vertex. 

(2) 

Near Ares , the numerator was small due to the proximity 
to the bottom of the band, and effectively reduced the 
divergence of the Raman vertex to <5~ 1//2 instead of <5~ 3 / 2 . 
The lower resonance frequency, Afil, however, is rather 
far from the bottom of the band so that the numerator 
does not possess any smallness. As a result, the strong 
enhancement of the Raman vertex as one approaches Afil 
from below turns out to be comparable, and for some val- 
ues of parameters even larger than the intensity near the 
high-frequency resonance. We will explicitly demo nstrate 



this feature in our numerical results in Sec. IV B 



C. Raman intensity at finite t' 

We now consider how the inclusion of a next-nearest 
neighbor hopping modifies the resonant behavior of the 
Raman vertex. We already discussed above that a 
nonzero t' breaks the particle-hole symmetry. In this 
case, the expression for the Raman vertex is more com- 
plex than Eq. ([?]) and has the form 



fdlk 



dk 



l^q^k—q ~ XqCk 



oj, - 2E k + 
1 



E 



k—q 



iT 



El + iT 



(30) 



where et and Ek arc defined as before and 
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E 



l k = -Atv k - 4t'cos(k x )cos(k y ), 

-3 — At' cos{k x )cos{k y ) . 



±s A 



the deviations from the resonance values, we obtain for 
G a (q) 



Here Eu' describes the energy dispersion of the conduc- 
tion and valence bands, respectively. 

One can easily see that at finite t', only two out of the 
three terms in the denominator in Eq.(31) can vanish si- 
multaneously. A nonzero t' thus effectively transforms 
the triple resonance into a set of double resonances. Ob- 
viously, there exist five combinations of terms for which 
the denominator can vanish. One of these, namely the 
one with Ui — 2E k = and ujf — 2E k - q — 0, yields a res- 
onance in exactly the same region of the (u>i, Alu) plane 
where the resonance occurs without t' . We will show that 
of the remaining four combinations only two are truly di- 

(2) 

vergent in the vicinity of uir es - 

In order to obtain some analytical results we again have 
to calculate the effective Raman vertex which now has a 
more complex form: 



G a {q) = 2vr y dx 



fde k \{Oe k - q u . 



de 



k q 
1 



(ui - 2E k + iT^j (u f - 2E k _ q + iT 



(Aw - 2uj q + iT 



Vi-v 9 + E%-E%_ q + ir 



El 



rpv 

b k-q 



iT 



We first observe that while the energy dispersion of the 
quasiparticles depends linearly on t' , the magnon disper- 
sion and the magnon-magnon scattering vertex V(k, q) 
depend only on (t') 2 . Since (t' jt) 2 < 0.25 (otherwise, the 
antiferromagnetic state is unstable), we will just neglect 
t' in uj q and V(k, q). 



Consider first the situation near 



,(2) 



when there is 



a true resonance in M R . Performing now the same 
manipulations as before, i.e., expanding near the bottom 
of the band and neglecting the numerator in Eq. (pl|) , we 



obtain M e R ff = (M e R n (a)+M e R n (-o))/2, where a = t'/J 
and 



M R ff (a) oc 



({ 2 i 



G a {q) 



(AA-gV2 + ir) 



(32) 



with 



G a (q) = J d 2 k (A, - 4 + 2k 2 + iT)- 1 

x (A.j - 4 - AA + 2{k - q) 2 + iT)- 1 
x 2X t — 8 - AA + 2fc 2 (l + a) + 2{k-qf{\ - a) 



+ (A\-qV2)+iT 



(33) 



Expanding near the point where the first two terms in 
the denominator in G a (q) vanish and integrating over 



1 



(x + iT)(k x(l + a) /{I - a) + z + iT) 
1 

. \/2k qCf cos 0o ~ (koq) 2 sin 2 (j) + iT 

^ • (34) 

J 2koqC% cos 0o — (koq) 2 sin 2 </>o + iT. 

Here k = [(4 - X l )/2} 1 / 2 and cos^o = (2q 2 - AA)/4fc g 
are the same as for the resonance with t' = 0, z = 
(AA(1 + a) - g\/2)/4(l - a), and Cf 2 are given by 

CI = x(k - qcos<f) ) ; 

C| = x(2k /(l + a) - q cos o z) . 

One can easily verify that the derivatives of Cf 2 are neg- 
ative in which case the poles and branch cuts in Eq.(p34) 
are located in different half-planes. This implies that the 
integral over x is finite. Performing the explicit integra- 
tion over x, we obtain after some simple manipulations 

Ga(q) = t— x 

z k q 



(31) 



sin 



90 



2z cos 0Q , -p 
k q 



1 



2z cos <pQ I 
k q 



I sin <p | 



(35) 



Now we are left with the integral over q in Eq.(|32|). In 
Sec. Ill A, the q integration was confined to a narrow re- 



gion around (ft = AX/y/2 and yielded Mjj cx \5\ 3 / 2 

where 5 — Xi — Arel- At finite t', an analysis of Eq.(|35|) 
shows that there exist two regions in momentum space 
which yield singular contributions to M R (a). The first 
region is still the vicinity of q = q a = AA/V2. How- 
ever, since z is finite for q = qo, this region yields a 
weaker, <5 -1 / 2 singularity in M R ^(a). In practice, ev- 
erywhere except for the immediate vicinity of the reso- 
nance, this divergence is fully compensated by the numer- 
ator in M R (a) which vanishes linearly as A^ approaches 

(2) 

ymax _ 4 which j s very close to Ares- We checked that 
the same functional behavior also holds for M R ^ (—a). 

The second singular contribution to M R ^(a) comes 
from the q integration over the region where z is nearly 
zero. The conditions z — and sin 4>o = specify a line 
in the (A,-, AA) plane with 



A a = A( 2 e )( a ,AA) = 4- (1 -y( 1 

4(1 + i 



a) 2 ) 2 



(36) 



We found that near this line, the Raman vertex also di- 
verges as <5 -1 / 2 where S now measures the deviation from 
Arei (a). This square root divergence near z — also 
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holds for Mftf(—a) for which the resonance incident fre- 
quency is given by Eq.(|36|) with a replaced by —a. We see 
therefore that a nonzero t' splits the strong resonance at 
Arel with a (5 -3 / 2 singularity into three weaker resonances 
with S^ 1 / 2 singularities. One of these weaker resonances 

(2) 

still occurs at A r es while the two new resonances occur at 

(2) 

A res (±a). For a = 0, the three reson ance lines coincide 



and we recover the result of Sec |III A 
For a 



-0.5, which is relevant to the cuprates, and 



(2) 

1.4 we have A res 



AA 

(2) 

Ares (—a) ~ 3.58. We see that the new resonance frequen- 
cies are further away from \™ lax than Aril and therefore 
should be less effected by the smallness of the numerator 
in Mjj . Naively, this should make the new resonances 

(2) 

stronger than the one at X res - However, we found that 

(2) 

the overall numerical factor is larg er near Ares than near 
the two new resonance lines. In this situation, t' just 

(2) 

reduces and broadens the peak at A re s without actually 
producing comparable peaks at the two new re sonance 
frequencies. Our numerical findings in Sec. IV C| are fully 
consistent with this result. 

Finally, we shortly discuss the effect of t' on the low- 
frequency resonance. We found that the inclusion of t' 
shifts the frequency range for the enhancement due to 
the branch cut but does not introduce any new physics 
near Aril. We again obtained that there is no real diver- 
gence of the Raman vertex in this frequency range but 
that slightly below Aril , the vertex acquires a branch cut 
enhancement which mimics the resonance behavior. The 
calculations near A r ll are, however, rather involved, and 
we did not succeed in fully solving the problem analyti- 
cally. We will discuss our numerical results near A r ll in 



Sec. IVC 



IV. NUMERICAL RESULTS 

In the following subsections we will present our nu- 
merical result s for the R aman lineshape and the TMPH. 
In Sees. IV A and IVB we first consider a system with 



particle-hole symmetry. In Sec. IV C we study how the 
form of the TMPH is modified due to a finite next-nearest 
hopping term t' whic h breaks the particle-hole symmetry. 
Finally, in Sec. IV D we discuss how the renormalization 



of the interaction between light and quasiparticles due 
to vertex corrections affects the form of the TMPH. We 
summarize all relevant formulae for the numerical compu- 
tation of the Raman intensity with final state interaction 
in Appendix 

Before we proceed, we want to point out the differ- 
ences in our numerical and analytical considerations for 
uji w u>res- For our numerical calculations we use the 
mean-field form of the fermionic excitation spectrum, 
which in the case t' — is degenerate along the boundary 
of the magnetic Brillouin zone. This particular form of 



the dispersion yields, besides an enhancement due to a 
branch cut, also a real divergence of Mr at will, though 
with a small overall factor. In our analytical calcula- 



tions in Sec. [II B we replaced this mean-field form by 
a quadratic dispersion around the top of the band, in 
which case the divergence transforms into a strong en- 
hancement. We therefore expect that our numerical re- 
sults will overestimate the strength of the low-frequency 
resonance. 

Finally, we shortly discuss some technical a spe cts of 
our numerical calculations. It follows from Eq.(Al) that 
the expression for the Raman intensity contains four- 
dimensional integrals with strong singularities. In order 
to make a numerical evaluation possible, one has to in- 
clude a fermionic damping, which cuts the singularities. 
However, if the damping is too large, subleading terms 
become stronger than the triple resonance effect. We 
found, for example, that a fermionic damping T = 0.4J, 
which was used in Ref. [E5| almost destroys the reso- 

(2) 

nance at uif-es ■ We therefore only consider relatively small 
fermionic dampings with 0.05J < T < 0.10J. Further- 
more, in order to ensure sufficient accuracy of the re- 
sults, we evaluated the necessary integrals on lattices up 
to 1000 x 1000 sites. We verified in each case that the 
convergence of the results was satisfactory. 



A. Raman line shape in B± g geometry 

We first present our numerical results for the Raman 
lineshape as a function of Aw for fixed w^. Our main 
result is that the Raman lineshape evolves with increas- 



es 



(i) 



to 



< W r el, an d 
LUrea- in 



ing u>i from a slightly asymmetric form at Wj 
a strongly asymmetric form at will < 
then back to an almost symmetric form at w 
order to show this we present the results for three inci- 
dent frequencies: w, ss 2A + 2.9J, « 2A + 6J, and 
u>i w 2A + 7.9J. In the first and third case the triple 
resonance and the two-magnon peak positions coincide, 
whereas in the second case they are well separated. 
) ~ ,W 

aj LOi ~ LiVes 

The Raman intensity as a function of transferred fre- 
quency without and with a final state interaction is pre- 
sented in Fig. |^a,b, respectively. The main difference 
between the two figures is the presence of an unphysical 
singularity in Fig. @a at Auj max = AJ which is due to a 
divergent density of states at the boundary of the mag- 
netic Brillouin zone. As in the LF-theory, the inclusion 
of a magnon-magnon interaction eliminates this singular- 
ity as is seen in Fig. ^o. A more relevant point is that 
both figures contain a strong peak at Aw — 2.8J. While 
the peak in Fig. ^|a is solely due to the divergence in the 
Raman vertex at will , the peak in Fig. ||b is a combined 
effect of the resonance in the Raman vertex and multi- 
ple magnon-magnon scattering. We see that the peak in 
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FIG. 6. The Raman intensity as a function of transferred 
OJrea a ) without and b) with a final 



frequency Aw for Wi 
state interaction. 



Fig. [|a is strongly enhanced by the final state interaction. 

Furthermore, we see that the Raman lineshape in 
Fig. ^|b is slightly asymmetric with a larger intensity at 
higher transferred frequencies. This asymmetry is most 
likely to be a property of the Raman vertex since the fi- 
nal state interaction yields a symmetric peak. One can 
indeed see this asymmetry already in Fig. ^|a. The two- 
magnon lineshape obtained numerically is consistent with 
our analytical results in Sec. [II B. There we attributed 



the asymmetry of the two-magnon profile to the branch 
cut in the Raman vertex which for u>i = 2A + 2.9 J exists 
only for Aw > 2.8 J. 

b) LOres < Wj < LOres 

The form of the Raman profile changes quite strongly as 
one moves from Wrel to intermediate incident frequencies. 
In Fig. [^we present, as an example, the Raman intensity 
including a final state interaction for u>{ — 2A + 6.0 J. A 
comparison with Figs. |^ shows that the anisotropy of the 
intensity is now much stronger. This result is quite ex- 
pected since in this range of incident frequencies, the Ra- 
man vertex resonates at transferred frequencies above the 
two-magnon peak. In particular, for w, = 2A + 6J, the 
triple resonance occurs near the maximum transferred 
frequency Aw = 4 J (see Fig. ||). 
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FIG. 7. The Raman intensity as a function of transferred 
frequency Aw for Wj — 2A = 6.0J in the interacting case. 



The two-magnon profile for intermediate frequencies 
within the triple resonance theory was earlier obtained 
by Schonfeld et ai. Our results arc in full agreement 
with theirs. 

) ~ ( 2 ) 
C) W l W LOres 

The results for the intensity with and without final state 
interaction are presented in Fig. ||b,a, respectively. The 
intensity without a final state interaction again exhibits 
an unphysical divergence at the maximum transferred 
frequency Alo = 4 J which disappears when one includes 

(2) 

a magnon-magnon scattering. Near uj r es, however, this 
divergence is confined to a very narrow region near 4 J. 
Furthermore, we obtain that in both figures the peak at 
around Aw = 2.8 J is almost symmetric. This is also con- 
sistent with our analytical results in Sec. Ill A. 
In addition to the peak at Aw = 2.8J, both intensities 



also possess a slight maximum around 3.3J which prob- 
ably originates from subleading, branch cut terms in the 
intensity. 

The evolution of the Raman profile with increasing 
u>i from slightly asymmetric form around w^' to a pro- 
nounced shoulder-like behavior for intermediate frequen- 

(2) 

cies, to a symmetric form close to LOres IS fully consis- 
tent with the experimental results on SriCuOiGl2 and 
Y BaiCu^O^.i. We consider this agreement with the data 
as yet another evidence that the triple resonance diagram 
dominates the scattering process in the resonance regime. 



B. Two-magnon peak height 

We now discuss the TMPH as a function of w,. For 
the calculation of the TMPH, we fix the transferred fre- 
quency at a value which corresponds to the maximum of 
the two magnon profile (which, depending on w^, occurs 
between Aw = 2.8J and Aw = 2.9J) and plot the inten- 
sity of the maximum as a function of w^. We present the 
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FIG. 



The Raman intensity as a function of transferred 
frequency Auj for uji — luIVs, a) without and b) with a final 
state interaction. 



results in Fig. ^ for two different values of the fermionic 
damping T. In both cases we clearly observe two max- 
ima at Wres ~ 2.9J and u^el ~ 7.9J. The positions of 
these maxima are in good agreement with the analyti- 
cal predictions and the experimental data. The form of 
the TMPH near uj^el is clearly asymmetric: the inten- 
sity drops faster above the peak than below. This form 
agrees with our analytical results. For intermediate inci- 
dent frequencies (4.0J < (w t - 2A) < 7.5 J) the TMPH 
remains basically constant and, in addition, is practically 
T- independent. This behavior, we believe, results from 
the fact that in this frequency range the triple resonance 
occurs at Aw = 4 J which is too far away from the two- 
magnon peak to influence its height. 

(2) 

Upon increasing T, we find that the TMPH around uir es 
drops much more rapidly than around Wrel ■ This is fully 
consistent with our analytical result that the divergence 
in the Raman vertex, which is only cut by the fermionic 

(2) 

damping, is present only near ujr es while the maximum 
near u4e S is just an enhancement which does not crucially 
depend on the damping. 

We also obtained two results which are not fully con- 
sistent with the experimental data. The first one is 
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FIG. 9. The TMPH as a function of incident frequency uji 



the ratio of intensities at the two maxima. We found 

(2) 

that while the divergence in Mr exists only near Wr es , 
the non-divergent terms are much stronger around uires- 

As a result, the ratio of intensities I (ivill) / 1 (luIVs) for 
r = 0.05J is sa 1, while experimentally, this ratio is 
clearly smaller than one, though the actual number dif- 
fers between Sr-iCuOzCli and Y BaiCu-sO§.\. The sec- 
ond discrepancy concerns the behavior of the TMPH in 

(2) 

the vicinity of u>r es . Analytically, we found that the 
TMPH should follow an inverse linear behavior at some 

(2) 

distance from ujr es and an inverse cubic behavior in the 

(2) 

immediate vicinity of Wrei- To verify this result, we cal- 
culated the TMPH for several frequencies in the vicinity 



of, 



,( 2 ) 



and present the results in Fig. nty (the dashed line 



in this figure is a guide to the eye). Within our numerical 
accuracy, we indeed found an inverse linear dependence 

(2) 

which, however, only exists for a small region near u>r e ' s , 

(2) 

namely for 0.1J < u))- es — uji < 0.25J. Experimentally, 
this region extends over a much wider frequency range of 

(2) 

about lev. Very close to ujr es , the divergence is cut by 
the fermionic damping, and it is impossible to verify the 
predicted inverse cubic behavior. 

We already mentioned that one of the reasons for the in- 
correct ratio of intensities I (u>res) / 1 \^?es) lies in the over- 
simplified mean- field form of the fermionic dispersion, 
and, in particular, in the degeneracy along the boundary 
of the magnetic Brillouin zone. One would thus expect 
a better agreement with experiments if this degeneracy 
is lifted e.g., by the introduction of a finite next-nearest 
hopping t' . We address this issue in the next subsection. 



C. Raman intensity for a nonzero t' 



In Sec. [II C we found that a nonzero t' splits the triple 

(2) 

resonance around uj res into three double resonances one 
of which occur in the same region of the (cj; 7 Auj) plane 
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FIG. 10. The inverse Raman intensity for Aui = 2.9J as a 
function of ajfes — tJt- 
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FIG. 11. The shaded area describes the region of the 
(wi,Aw) plane in which a double resonance occurs for 
u, - 2E k = and Wi - ^ + - = 0. The solid 

line corresponds o q x — q y and the dashed line to q x — 0. 



as the triple resonance in the absence of t' while the other 
two occurs in different regions of the (wj, Alu) plane. In 
Fig. [ll] we plot the region of the (u>i, Auj) plane in which 
one of the remaining double resonances occurs. The form 
of the shaded area is similar to Fig. 0. We see that in 

(2) 

the vicinity of cOres the new resonant region reduces to 
a single line, just as the resonance for t' = 0. Near w£es 
the situation is more complex since the different double 
resonances overlap. 

In Fig. ||we present the result for the TMPH for t'/t = 
-0.16. A comparison with Fig. gfor the TMPH at t' = 
shows that a finite t' reduces the TMPH at both reso- 

('12). 

nant frequencies tore's ■ This reduction is fully consistent 
with our analytical calculations since now only two terms 
in the denominator of the Raman vertex vanish simul- 
taneously while the third scales as 0[t' /t). The reduc- 
tion, however, is not uniform, and the TMPH around the 
high frequency resonance decreases much more rapidly 
than the one around the low-frequency resonance. Most 




0); - 2A (in units of J) 

FIG. 12. The TMPH as a function of the incident frequency 
uji for t'/J — —0.3 and T = 0.05 J. The inset shows that 
despite the strong reduction one can still observe a maximum 
around w^es- 



probably, the increase of the ratio is caused by two ef- 
fects. First, the regions of double resonance overlap 
around Wres, but not around Wrel- Second, a nonzero t' 
also affects the interaction vertex Vi / between light and 
fermions and reduces it much more strongly around ujr es 
than around uiiW. To see this, we recall that in the mean- 
field approximation we have Vif = (dek/dk) §ij with 

(2) 

ej; = — 2t (cos k x + cos k y ) + 4\t | cos k x cos k y . Near uires , 
the dominant contribution to the Raman vertex comes 
from fermions near the bottom of the band (k sa 0) in 
which case the vertex between light and fermions is re- 
duced by a factor of (1 — 2j£'j/t). In contrast, the reso- 
nance near u^X is dominated by fermions near the top of 
the valence band (k = (±7r/2, ±7r/2)) in which case the 
effect of t' is negligible. 

We see therefore that the inclusion of a finite t' actu- 
ally worsens the agreement with the experiments since 
the ratio of intensities increases. In the next subsection 
we will consider whether vertex corrections can possibly 
reverse the effects of t' and restore the correct quantita- 
tive behavior of the TMPH. 



D. Vertex Corrections 

There are several vertices in the diagram for the Ra- 
man matrix element, each of which is renormalized by 
vertex corrections which are generally not small at large 
U. The calculation of all vertex corrections is beyond our 
computational abilities and in this section we will there- 
fore focus on the corrections to the vertex between light 
and fermionic quasiparticles, Vif. Some evidence that 

(2) 

the vertex between light and fermions near LUr es is larger 
than in the mean-field theory comes from the measure- 



15 




Two-Magnon Peak Height 



FIG. 13. The lowest order vertex corrections to the 
light-quasiparticle interaction. The solid and dashed lines 
represent the conduction and valence band quasiparticles, re- 
spectively. The internal wavy line describes the exchange of 
a magnon. 
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ments of the optical conductivity in Gd,2CuOi, PriCuO^ 
and YBa^CujP^ pSfl . These experiments have demon- 
strated that the measured conductivity is larger than 
the one calculated with the mean-field form for Vi / even 
though it basically follows the same frequency depen- 
dence. We will study the vertex corrections to Vij semi- 
phenomenologically and our goal will be to illustrate how 
they can, in principle, reverse the effects of t' . 
The lowest order correction to V// in a formal perturba- 
tive expansion in 1/S is presented in Fig. [l3| A simple 
analysis shows that the relative vertex correction scales 
as U / JS, i.e., it is small only in the limit of a very large 
spin. For realistic S, however, we have U/ JS ^> 1, and 
the corrections to Vif are large. This clearly implies that 
one should sum up an infinite series of corrections to ob- 
tain the proper renormalization of the vertex between 
light and fermions. We will not do this but rather model 
the effect of the vertex corrections phenomenologically 
by introducing an effective vertex in the form 

yf (A) = (^,,)(l + a^) (37) 

with a as a parameter. The effective vertex in Eq.(|37]) 
still possesses the same symmetry as the bare vertex 
Vif(k) and therefore still vanishes at the bottom of the 
band. However, the slope of Vf/^ (k) around k = can 
now be quite different from the mean- field result. 

We computed the TMPH with the effective vertex us- 
ing various values of a. The result for a = 0.5 is pre- 
sented in Fig. [l4]. A comparison with Fig. || shows that 
the effect of the vertex correction is rather strong; the ra- 
tio /(u/fes)/ 'I{^res) is decreased by a factor of about 2.5. 
In addition, we also observe a relative increase of the 
TMPH for intermediate 4.5J < (ua - 2A) < 5.5J. This 
last effect leads to an extension of the region in which the 
Raman intensity possesses an inverse linear behavior. 

The decrease of the ratio of the intensities and the 
extension of the frequency range of the inverse linear be- 
havior are both in agreement with the experimental re- 
sults f2(i|j . We therefore see that by adjusting the vertex 



FIG. 14. The TMPH for V = 0.05J and a = 0.5. 

corrections to Vif(k) without violating the symmetry re- 
quirements of the model, one can, in principle, not only 
obtain good qualitative, but also quantitative agreement 
with the experimental data. The question is, however, 
whether, e.g., a = 0.5 which we used in Fig. can be 
obtained in a microscopic calculation. These studies are 
clearly called for. 

V. DISCUSSION 

We first summarize our results. The intent of this pa- 
per was to study the full Raman intensity in the reso- 
nant regime by simultaneously considering the effects of 
the triple resonance in the Raman vertex and the final 
state magnon-magnon interaction. We derived an ex- 
plicit expression for the full Raman intensity in the reso- 
nant regime as a function of both, transferred frequency 
Auj and incoming frequency Wj. We obtained analyti- 
cally and numerically the two-magnon Raman profile as 
a function of the transferred photon frequency Auj and 
the dependence of the two-magnon peak height on the 
incident photon frequency Wj. We found that the reso- 
nant behavior of the Raman vertex survives the inclu- 
sion of a magnon-magnon interaction and obtained two 

maxima in the peak height at coles ~ 2A + 2.9 J and at 

(2) 

aires ~ 2A+7.9 J. The position of the two maxima are the 
same as in the semi-phenomenological approach by CF 
which considered the triple resonance enhancement and 
final state interaction independent of each other. We first 
studied in detail the two-magnon profile at various inci- 
dent frequencies. We found that the two-magnon peak 
is slightly asymmetric near w^es with larger intensity at 
higher frequencies. As the incident frequency increases, 
the asymmetry becomes stronger, and the two-magnon 
profile acquires a shoulder-like feature above the peak. 

This is consistent with earlier results p5| . For frequen- 
ts) 

cies around u>re S , however, we found that the anisotropy 



1G 



disappears, and the Raman profile acquires almost the 
same form as in the nonresonance, LF regime. 

We then proceeded to a more detailed study of the 
two-magnon peak height. We verified that the inverse 
linear behavior of the Raman intensity near uir ea sur- 
vives the effect of the final state interaction. Further- 
more, we considered the behavior of the Raman vertex 
near cures ■ We found in our analytical considerations for 
which we assumed an isotropic dispersion near the top of 
the band that the divergence is almost completely sup- 
pressed. However, the Raman vertex contains a branch 
cut which gives rise to an enhancement of the intensity 
in some range of frequencies u>i < toiW which termi- 
nates only slightly below u>res ■ In our numerical calcula- 
tions, for which we considered a mean-field form of the 

dispersion, we obtained a weak singularity at cjrel but 
also a strong enhancement of the Raman intensity for 
u>i slightly smaller than u>le\- This last enhancement is 
virtually independent of the damping. 

We found that the ratio of the Raman intensities 
7(wres)/-f(wres) is already rather large for small damping, 
contrary to the assertion by CF. A much smaller ratio is 
needed for a quantitative agreement with the experimen- 
tal data. We attribute the large ratio to an unexpectedly 
strong enhancement of the two-magnon peak near uj^Js 
due to a branch-cut anomaly in the Raman vertex. 

We further studied how the the triple resonance is 

( 2) 

modified by a next-nearest hopping term t' . Around uj res , 
we found that the triple resonance is split into three dou- 
ble resonances, but the linear divergence of the Raman in- 

(2) 

tensity near uir es is not changed. This splitting, however, 

(2) 

reduces the intensity around u! res relative to the intensity 
around wi-el where the effect of a finite t' is rather weak. 
As a result, the ratio of the intensities I (uji]) s ) / 1 (ui^l) 
increases. 

Finally, we have demonstrated that the ratio of the 
intensities at the two resonance values of u>i is sensitive to 
the actual form of the vertex between light and fcrmions. 
We have shown that the corrections to the mean-field 
vertex are large and modeled their effect by introducing 
an extra factor (1 + az/|) into the vertex. We considered 
a as an adjustable parameter and showed that the ratio 
of the intensities can be substantially reduced already for 
moderate values of a. 

We now discuss our results in the context of the key 
experimental features that we listed in the introduction 
as being in disagreement with the LF theory: 
1) Changing lineshape with oji 

Our results for the evolution of the Raman profile with u>i 
is in complete agreement with the experimental results by 
Blumberg et al. [|fj on Y BaiCu^O^. For Sr 2 Cu0 2 Ch, 
the highest experimentally accessible frequency is smaller 

(2) 

than ujres and we therefore do not know whether the 

(2) 

Raman profile eventually becomes symmetric near ujr es . 
For intermediate Ui, however, our results agree with the 



experimental data. 

An issue which we have not addressed in our ap- 
proach is the actual rather than relative width of the 
two-magnon peak. Experimentally, it is much broader 
than in our model. Previous studies by Weber and Ford 
p7f , however, have shown that the broadening may be 
due to a magnon damping. They demonstrated that a 
small damping due to, e.g., an interaction with phonons 
already gives rise to a considerable broadening of the two- 
magnon peak. This result has also been obtained in nu- 
merical studies p8| . 

Another reason for the broadening of the Raman line- 
shape is the fermionic damping T. Incidentally this 
damping may also account for the experimentally ob- 

(2) 

served additional broadening of the TMPH around Ur es 
since in this frequency range, the dominant contribution 
to Mr comes from fermions with k w which exhibit the 
largest damping pOfl . 

2) The TMPH as a function of u>i 
The two key experimental results for the TMPH, we re- 
mind, are the presence of two maxima in the TMPH, of 
which the higher frequency maximum is stronger in all 
compounds, and an inverse linear behavior of the Raman 

(2) 

intensity near the upper resonance frequency ujr es . In 
our analytical and numerical calculations we found the 
two maxima in the TMPH whose positions fully agree 
with the experimental data. In addition we found that 
the low-frequency maximum in the TMPH is anisotropic 
with a higher intensity at lower frequencies which is also 
consistent with the experimental results. 

Our numerical data, however, differ quantitatively 
from the experimental results in that the ratio 

J(wres)//(wrel) is too large. On the basis of our an- 
alytical results we would expect the opposite behavior 
since we found that the actual resonance in the Raman 
intensity (i.e. a divergence in the absence of a fermionic 
damping) exists only near w^el while the peak near uiW 
is just the enhancement due to nonsingular terms in the 
Raman vertex. It turns out however that these nonsin- 
gular terms are anomalously large. Naively, one would 
expect that the inclusion of a next-nearest neighbor hop- 
ping t' would lead to an improved agreement of the ratio 
with the experimental data. In contrast, we found that 
a finite t' suppresses the high-frequency resonance even 
further. On the other hand, we have demonstrated that 
the inclusion of the corrections to the interaction between 
light and fermions may substantially increase the vertex 
near ll>1 2 J s compared to the vertex near Wrel- This eventu- 
ally yields a much better ratio of I (u>res) fibres) which 
can be made fully consistent with the experimental data 
by adjusting the magnitude of the vertex correction. 

Our analytical and numerical computations also re- 
produced the inverse linear behavior of the Raman in- 
tensity, which was observed in Y BaiCu-zO§.\ [Q and 
PrBa-iCuzO'j and, to a certain extent, also in 

Sr^CuOiCli. However, we also found that the range 
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FIG. 15. The Raman intensity as a function of transferred 
frequency Auj for constant denominator. 



of uji in which this behavior was observed experimentally 
is much larger than in our analysis. The inclusion of the 
vertex corrections improves the agreement with the data 
but does not make it perfect. This issue requires further 
study. 

In conclusion we have provided a detailed study of Ra- 
man scattering in the resonant regime. We confirmed 
that the key experimental features of magnetic Raman 
scattering can be explained qualitatively, and to some 
extent quantitatively within the triple resonance theory. 
We believe that the remaining quantitative discrepancies 
are due to insufficient knowledge of the quasi-particle en- 
ergy dispersion, lifetime effects and the form of the vertex 
function between light and fermions. 

A final remark. We discussed in Sec 
intermediate incident frequencies, 



UJres 



IV A 



that for 

<. Lu <. LiVes, 

the triple resonance occurs relatively far from the two- 
magnon peak and, to first approximation, does not in- 
fluence the two-magnon lineshape. In other words, in 
calculating the two-magnon lineshape, one can, with rea- 
sonable accuracy, set the denominator in the triple res- 
onance diagram to a constant and compute the Raman 
vertex in the same way as in the LF theory. It turns out 
that this procedure yields good agreement with the data 
not only in B\ g but also in other scattering geometries. 
To demonstrate this, we present in Fig. |l5| the results 
for the intensity in four different scattering geometries. 
These results have to be compared with the experimental 
data for SV^CuC^C 1 ^ and Y Ba^Cu^O^.x |3j|. which we 
reproduce in Fig. [l6| Since J ~ 1000cm , the compar- 
ison with Fig. |l5| is valid only for Auj < 4000cm -1 . A 
finite scattering intensity at larger Auj is probably due 
to multi-magnon scattering. We consider the agreement 
between the two figures as rather good and view it as an 
another piece of evidence in favor of the triple resonance 
theory. 
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FIG. 16. Two-magnon Raman scattering spectra from 
Sr2Cu02Cl2 and YBa2Cus06.i at room temperature in differ- 
ent scattering geometries. Courtesy of the authors of Ref. |}9| . 
The labels indicate: xx (B\ g + Ai g ), x'y' (Bi g + A2 g ), x'x' 
{Aig + B 2g ) and xy {B 2g + A 2g ) The excitation energy is 
2.73 eV, which is about 0.4 eV below uii-es- 
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APPENDIX A: THE RAMAN INTENSITY WITH 
MAGNON-MAGNON INTERACTION 

In this appendix we present the formulae for the nu- 
merical computation of the full Raman intensity. Our 
starting point is the Golden Rule formula, Eq.(Fy) 



,Aw)cx J 0|A4 ot (^,Ao;,< ? )| 2 ( 5(A^-2a; g ) 



(Al) 



where M^ ot is diagrammatically presented in Fig.^j. The 
Golden Rule formula for the intensity corresponds to the 
diagram in Fig. [l]a in which the intermediate magnons 
are on the mass shell. 

The analytical expression for M|j ot has the form 
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M R ot 



1 + 7/45 



(A2) 



ij) , and Mr and are given by 
and by Eqs.© and (|lj) for 



where 7 is given in Eq. 
Eqs.(g) and (§) for f : 
t' 7^ 0, respectively. 
Note that we use the full form of A7r and do not project 
it on v q as was done in Ref . p5fl . Our numerical computa- 
tions show that especially for small T, the B± g component 
of Mr has a more complex dependence on the magnon 
momentum than just v q . 
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